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SUMMARY 


This is the final report presenting the results of the analytical 
and experimental investigations of flow fields downstream of annular jets, 
both with axial and swirling flows. 

An analytical model was developed to predict the base pressures 
established by axial or swirling flow through annular Jets having blunt 
bases. The model made predictions for configurations having various flow 
angularities, radius ratios and swirl distributions. Approximately 50 
hours of scale model testing were conducted to determine several empirical 
constants that were required for the analytical model. 

The scale models tested consisted of 11 sets of interchangeable inner 
and outer exhaust flowpaths encompassing flow angles of 10 degrees radially 
outward to 20 degrees radially inward and radius ratios of 0.40 to 0.70. 

All of the models were designed to have a constant nozzle exit area of 
6.41 square inches (41.35 square centimeters). Various model configurations 
were tested with: 

• Axial flow 

• Swirling flow 

• Variable pressure ratio 

• Simulated louvers (blockage) 

• Skewed flow distribution 

Having evaluated the necessary empirical constants using the experi- 
mental results from the configurations tested with axial and swirling 
flows, a comparison of the theory with test results was made to verify the 
selection of the constants. A fair agreement was obtained considering 
the relatively small sampling of test configurations and the complexity 
of the analysis. Based on the comparison of test and analytical parameters 
such as velocity, flow and pressure coefficients, the analytical model 
appears to be an adequate representation of the exhaust flow downstream 
of annular nozzle with blunt bases. 



INTRODUCTION 


During the initial testing of a rotor-alone turbotip lift fan system, 
a deficiency in fan thrust was observed and identified as an excessively 
low hub base pressure. A program was initiated under NASA Contract NAS2- 
5462 to investigate, both analytically and experimentally, the effects of 
various types of flow fields upon the base pressures of an annular jet. 

This report covers this investigation. 

An analytical model was developed to predict the level of base 
pressures which would exist for various flow angularities and nozzle 
exit radius ratios of annular jets. The analytical model was supported 
by approximately 50 hours of scale model testing conducted as a General 
Electric Company, Evendale test facility. Empirical constants from the 
experimental results were evaluated and used to verify the analytical 
results with the experimental results. Two of the test configurations 
were used as an independent check on the validity of the empirical constant 
selection. 

Based on the comparison of test and analytical parameters such as 
velocity, flow and base pressure coefficients, the analytical model appears 
to be an adequate representation of the exhaust flow downstream of annular 
nozzles. 
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ANALYTICAL INVESTIGATIONS 


ANALYTICAL MODEL 

General Discussion 

When flow exits from an annular nozzle with a blunt centerbody, the 
pressure over the aft face of the centerbody may be considerably different 
from the ambient pressure into which the jet is exhausting. This base 
pressure may have a significant effect on the flow and thrust coef f icieVits 
of the nozzle. If the jet exits axially without swirl base pressure co- 
efficients, based on average jet velocity head, on the order of -.10 to 
-.15 are typical. However, if the jet has radial or circumferential velocity 
components, it is found that the base pressure coefficients may vary widely 
from these values. This analysis was undertaken to develop an analytical 
model which could be used to predict base pressure, thrust, and flow co- 
efficients for annular nozzles with nonaxial exit flow. The analysis also 
considers the effects of radial distribution of swirl, nozzle radius ratio 
and nozzle pressure ratio. 

Conceptually, the reduction of the base pressure below ambient pressure 
can be thought of as being the result of three separate factors. First, 
in a jet with axial exit velocity, the shear forces between the high velocity 
main stream flow and the relatively stagnant air mass composing the center- 
body wake must be balanced by an equal and opposite force acting across the 
base area. This component of base pressure is a true drag force and results 
in a thrust loss. 

Second, there is a component of base pressure due to the meridonal 
curvature of the flow streamlines in the vicinity of the nozzle exit. This 
streamline curvature causes a radial static pressure gradient in the flow 
and, thus causes the base pressure to differ from ambient pressure. If 
the flow exits axially, the streamline curvatures are small; but if the 
flow is angled radially inward or outward, the streamline curvatures are 
much larger, and the effect on the base pressure may be significant. This 
component of the base pressure has no effect on the nozzle velocity co- 
efficient since it is essentially a potential flow phenomenon. 
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Thirds if the Jet has a swirl velocity component, a radial pressure 
gradient must exist in the flow stream to support the centrifugal forces 
generated by the swirl. This then results in a base pressure lower than 
the ambient pressure. Furthermore, the swirling jet induces a rotational 
motion to the centerbody wake which in turn causes a radial pressure 
gradient across the base area and a further reduction in average base 
pressure. This component of base pressure can cause a large loss in nozzle 
thrust coefficient and represents the energy lost in the swirl velocity 
component. 

It is evident that the above three effects are closely interrelated so 
that they cannot be treated independently in a mathematical analysis of 
the problem; however, it is conceptually enlightening to recognize the 
separate factors affecting the base pressure. 
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ANALYSIS 


Assumptions and Boundary Conditions 

Figure 1 shows the general scheme of the flow model assumed in this 
analysis. The flow downstream of the nozzle exit is divided into two 
distinct regions, the main stream flow and the centerbody wake^ The 
dividing line between these two regions is the hub streamline of the main 
flow. It is recognized that there is actually flow interchange across 
this boundary, but it can be taken as the line across which the time 
averaged flow is zero. Thus, the total time averaged mass of air within 
the wake region is constant. It is assumed that static pressures In the 
two regions are equal along this streamline and that shear forces between 
the two regions may exist along this line. 

The outer streamline of the main flow is assumed to be an isobaric 
surface where the static pressure equals the ambient pressure into which 
the jet is exhausting. Shear stresses and mixing along this outer boundary 
are assumed to have negligible effect on the base pressure of the centerbody. 

At Station 1, the upstream boundary of the main flow, it is assumed 
that the radial distributions of total pressure, total temperature and 
angular momentum can be specified as will be explained later. At Station 
2, the downstream boundary of the main flow, it is assumed that the slope 
and curvature of all streamlines are zero. Total temperature and swirl 
angle distributions at this station are derived from the upstream values 
assuming that energy and angular momentum are conserved along each stream- 
line. The total pressure distribution at Station 2 is derived from the 
upstream total pressure and the assumption that the loss due to the shear 
stress between the main flow and wake is distributed uniformly across the 
main stream. 

Other assumptions are that the flow is steady, axisymetric and has the 
properties of an ideal gas. At the inlet station, it is assumed that the 
radial distribution of angular momentum can be sufficiently well represented 
as a linear relationship to the stream function. 


- 5 - 



General Approach 


The general approach used in this analysis is as follows: 

First, it is assumed that the shape of the 50% streamline can be 
represented by an exponential equation of the form. 

R„, = [A+B (Z/L) + C (Z/L)^] e + R^2 (1) 

readily seen that A and B are determined by the initial radius 
and slope of the 50% streamline. The constants C and are determined 
empirically and will be explained in more detail later. L is a length 
parameter that controls the rate at which approaches its asymptotic 
value, Rjj^ 2, and will be adjusted during the solution to satisfy axial 
momentum requirements. 

• Starting with an initial guess for L, the radius, slope and curvature 

of 50/4 streamline are calculated at a series of axial locations starting 
with Z = 0 and continuing until approaches R ^2 sufficiently close. 

• An Initial guess for the mass flow rate, M , is then made, and at each 
of the above axial locations, a calculation is made along a line normal 
to the 50% streamline to determine flow properties and radial locations 
of the other main flow streamlines. Appendix I discusses the equations 
used for this calculation and the computer program section of this report 
discusses some of the calculation details. Briefly, however, good 
approximations to the streamline curvatures are made from the mean 
streamline slope and curvature and then the equations of momentum, con- 
tinuity and energy are applied in a direction normal to the flow while 
satisfying the boundary condition of ambient static pressure on the 
outer streamline. These calculations begin at the nozzle exit and 
progress downstream and add an incremental total pressure loss at each 
step to account for the centerbody wake mixing loss. This incremental 
pressure loss is calculated from the local hub streamline velocity, 
density and incremental mixing surface area. This loss is assumed to 

be distributed uniformly through the main flow at each axial station 
however . 
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• Next, using the static pressures and velocities calculated above, 
axial components of the pressure and momentum forces are integrated 
across the upstream Station 1 and the downstream Station 2, The 
pressure forces across the wake area are calculated using the hub 
streamline static pressures and an empirical correction based on the 
level of swirl at the hub streamline. The momentum forces across the 
wake areas are zero since the net amount of fluid in the wake does not 
change. 

• A check is then made on the overall balance between axial pressure and 
momentum forces acting on the fluid contained between Stations 1 and 2. 
The radius ratio at the nozzle exit station is also compared to the 
desired radius ratio. If the axial force balance and radius ratio are 
not close enough, new estimates are made for L and M and the procedure 
is repeated until the parameters are in agreement within a specified 
tolerance. 

It is readily seen that the mass flow rate, M , has a first order effect 
on the radius ratio. In general, the length parameter L has its first 
order effect on the force balance since it directly affects the mean 
streamline curvature and the surface area of the wake-main flow mixing 
surface* Both parameters do, however, have significant influence on 
both radius ratio and the force balance and the iteration scheme must 
recognize this. 

• After values for L and M are found which simultaneously satisfy the 
desired radius ratio and the axial force balance, then the desired 
values for flow coefficient, thrust coefficient, base pressure coeffi- 
cients, etc., are calculated from the converged solution. 

Mean Streamline Shape 

An exponential equation in the form of Equation 1 was selected to 
represent the shape of the 50% streamline. This equation can be made 
to have any initial radius and slope at Z = 0 and asymptotically approaches 
zero slope and curvature far downstream as we would expect of the real flow. 
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The length parameter L controls the rate at which the mean streamline 
radius approaches its asymptotic value of R^2 in effect, then controls 

the rate at which the centerbody wake washes out. It can easily be shown 
that A and B are determined by the desired initial radius and slope and are 
given by 

A = 1 - Rm2 (2) 

B = A+Ltan0^j^ (3) 

if we assume an initial radius, of 1. 

The third terra in Equation 1 whose raagnitude is controlled by the 
constant C is not needed to raatch the end conditions but was added to 
achieve a more realistic axial distribution of curvature along the mean 
streamline. If C is set equal to zero, the third derivitive of from 
Equation 1 is very high when Z is small. This implies a high axial 
gradient * of mean streamline curvature and thus a high axial gradient of 
static pressure in the centerbody wake. Initially, it was felt that 
selecting C so that the third derivitive of R^ would be zero at Z = 0 
would result in a more realistic curvature distribution for the mean 
streamline. To achieve this then it can be shown that 

C = (-A + 3B)/6 (4) 

Finally, in the comparing the calculated mean streamline shapes with 
the model test data, it was found that a better correlation could be 
obtained if C was defined as 

C = (-A + 3B + Cg)/6 (5) 

where was detennined experimentally and discussed in the Application of 
Analytical Methods Section of this report. Cg then controls the initial 

value of the axial curvature gradient of the mean streamline and has the 
same sign. 

The value of R^gi is indirectly based on experimental data. The 
Application of Analytical Methods Section of this report describes how the 
inner streamline radius at the downstream station was found to vary in the 
test program. If the inner streamline radius, is known at the downstream 
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station^ it is a relatively simple matter to calculate the mean streamline 
radius knowing the flow, swirl distribution, etc, and assuming zero slope 
and curvature at all streamlines. 

Station Calculations Normal to Mean Streamline 

Streamline Curvatures 

In order to apply the momentum equation in the direction normal to 
streamlines, the streamline curvatures must be known along a potential 
line. Since we are assuming a shape for the mean streamline, we know only 
its slope and curvature and must estimate a curvature for the other stream- 
lines. If we assume that the flow area between the mean streamline and any 
other streamline is approximately constant over a short distance, then we 
may superimpose the two cases illustrated in Figure 2. This then will 
approximately relate the curvature of any streamline to the mean streamline 
slope and curvature on the same calculation line by the following equation: 


l/r = 1/ [r„ + (R„ - R)/cos Bn, ] 


+ (R 


stn 


C6) 


Summary of Other Equations 

In addition to the above equation for streamline curvature, five more 
equations are derived in Appendix I which are required to define the main 
stream flow properties along a calculation line normal to the mean stream- 
line. These additional five equations involve the six unknowns - r. RV . 
p , T, R and \{n and are summarized below: 


RVq = (RVq )^ + X M (1 -Y ) 

P/Pn,^ = (P/P^q)(T/T^q) 


(II) 

'TO ^ ^ 

t/Tto “ 1 + IQ (is) 


^ ^ M cos 


(16) 


2tt p R 


B V, 


m „ 


M 


3Y 


2tt p R r 


+ MJ. r YtI. Q 


'm 


2 (RVg).^ 


(RVq) ^ (110) 
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Since the last two of the above equations are partial differential 
equations, they must be integrated along the calculation line. This was 
done using a second order Runge Kutta procedure and is described in more 
detail in the Computer Program Section of this report. The other four 
above equations are used as auxiliary equations in the integration 
procedure. 

Axial Force Balance 

If we remember that the outer surface of the main stream is at ambient 
pressure and that there is no net flow into or out of the wake region, we 
may write the following equation expressing equilibrium between the 
resultant axial pressure force and the change in axial momentum of the main 
stream. 

. } %1 Rt2 

“ d'f = 2 tt R(Pa - dR + 2rr jf R(Pa-P2^ 

The two terms on the right-hand side of Equation 7 may each be written 
as two terms, one of which represents the pressure force in the wake region 
and the other of which represents the pressure force in the main stream 
flow as follows: 

M dR + 2 tt ^ R(P^-P^) dR 

^^ H 2 -%2 

+ 2 rr y R (Pa-P 2 ) dR + 2 fr y R (Pg-P 2 ) dR 

0 Rjj2 


In this analysis^ the first Integral on the left is calculated as 
follows: ^ ^ 

“ -i ^'^zr’^z2> ‘^’1' = M ^ cos (9) 

Where it has been assumed that 


and 


B ^ B 
®2 = 
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( 10 ) 


the first integral on the right is calculated by 

R(Pa-Pi)dR = ttR^i ^2 Pl^)' ] 

HI 


where is an empirical coefficient which accounts for the effect of 
sWirl in the wake region on the hub base pressure. The third integral 
on the right of Equation 8 is calculated by 

R(Pa-P2) d* - n 1^,2 (P. - P^j) (U) 


Collecting Equations 8, 9, 10 and 11 and rearranging gives the final 
axial force balance equation used in this analysis; 


^ / ('ml ®nr''„2> ^ R(Pa-Pi> 

^ 0 \l 


dR 


:t2 


. f R(Pa“I’2) dR - ^ 

«'%2 

■ ^ *h 2 (Pa-PH2> = 0 




2 

i 

HI 


( 12 ) 


Overall Performance Calculat ions 


After a converged solution to the flow problem is obtained^ we must 
still calculate flow^ thrust^ velocity and base pressure coefficients 
since they are the nozzle characteristics of primary interest. To do 
this ^ the following procedure was used: 


• Firsts a mass averaged ideal jet velocity, was calculated using 
total pressure and temperature at Station 1 expanding to ambient 
pressure. Also, a corresponding ideal density, p^, was calculated. 

# An ideal mass flow was calculated as 


M. = p. A 
i n i 

the geometric 

angularity factor applied. 


(13) 

where A^ is the geometric annular area at the nozzle exit with no 


• An ideal thrust was calculated as 


F. 

1 


= «i 


V. 

1 


(14) 
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*"P2 ~ 


- V 


( 19 ) 


Thus, Cp^ is the base pressure coefficient at the outer radius of 
the centerbody and Cp2 represents the average base pressure co- 
efficient over the centerbody area. Both are based on an average 
nozzle exit dynamic pressure. 
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COMPUTER PROGRAM 


General Description 

The computer program was written to carry out the iterative solution 
to the equations presented in the previous section. Time-sharing Fortran 
was used. In general^ input to the program consists of: 

• Geometric parameters describing the nozzle exit, 

• Total pressure, total temperature, and swirl distributions at the 
nozzle exit. 

• A wake shear stress factor. 

• Three empirical constants. 

• An initial guess for L. 

• Parameters used to calculate rotor loading coefficients if the swirl 
is generated by a rotor. 

• Two control indicators. 

All input is supplied at the time-sharing terminal at the time of 
program execution. Output of the program consists of; 

• A general description of the flow field. 

• Nozzle base pressure, flow, velocity and thrust coefficients. 

• An approximate description of the rotor blade environment and loading 
coefficients if the swirl is generated by a rotor. 

Figure 3 is a schematic block diagram of the program showing the 
general calculation procedure and flow of information. The solid lines 
show information flow between main program elements while the dotted 
lines show information flow between the main program and subroutines or 
between subroutines. Numbers in the upper left-hand corner of some 
blocks refer to statement numbers in the program listing given in 
Appendix III, 

There are two options in running this program - the rotor option 
and the stator option. For the stator option, it is assumed that the 
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swirl is generated by a "no loss" stator with uniform inlet total 
pressure and temperature. For the rotor option, it is assumed that the 
swirl is generated by a rotor with uniform upstream total pressure and 
temperature and zero upstream swirl. The energy addition process is 
assumed to be isentropic but will result in a radial variation of total 
pressure and temperature at the nozzle exit station if a nonconstant 
angular momentum is specified at this station. 

"ROTFLO" Subroutine 

A key element of this program is the subroutine "ROTFLO" which 
performs the station calculations normal to the mean streamline. This 
subroutine uses a second order Runge Kutta method to simultaneously 
integrate Equations 16 and 110 along a station calculation line. Equa- 
tions II, 14 and 15 are used as auxiliary equations during this 
integration. 

"ROTFLO" can either begin with the radius given for the inner 
streamline and integrate outward to the tip streamline or begin with 
the radius given for the mean streamline and integrate both ways to the 
tip and hub streamlines. In either case, the initial value of the 
merldonal velocity, is Iterated upon until the static pressure on 
the tip streamline equals the ambient pressure. Twelve streamlines 
are used in the integration procedure and are, located at 0, 4, 10, 20, 
30, 40, 50, 60, 70, 80, 90 and 100 percent of the flow starting at the 
hub streamline. The streamlines are spaced closer together in the inner 
portion of the flow since the flow tends to vary more rapidly in this 
region. 

Iteration Scheme 

In order to arrive at a solution, both the mass flow rate, M, and 
the length parameter, L, must be iterated upon until the nozzle exit 
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radius ratio and axial force balance equation are satisfied. The scheme 
selected to do this has been called the secant method and will be briefly 
described for the case of two variables. 

Let be the percent error in radius ratio and be the percent 

error in the axial force balance and plot these errors against M and L 
for three previous iterations as shown in Figures 4a and 4b, These 
points lie on error function surfaces. If we pass a plane through each 
of these two sets of points, we will have ''secant” planes to the two 

error function surfaces. Suppose that in Figure 4a the secant plane 

• • 

intersects the M-L plane in the Line ab. Combinations of M and L which 
lie on Line ab can usually be expected to reduce the error ej^ to a lower 
value than the previous tries. Similarly points which lie along line 
cd in Figure 4b can be expected to reduce the error 62 # A combination 
of M and L then which lies on both ab and cd should then be an improved 
estimate for the next iteration. This then is the scheme which the 
Subroutine "INT2” uses to give improved estimates for M and L at each 
step in the Iteration. 

m 

In order to start the iteration procedure, initial estimates of M 
and L are made and then each in turn is perturbated by a small amount 
in order to give three starting points. 

Input Definitions 

Nozzle exit radius ratio, dimensionless. 

Mean streamline angle from axial; deg. + = outward, -- = 
inward. 

Wake mixing shear stress coefficient; based on hub streamline 
meridonal velocity head; dimensionless. 

Initial guess for the length parameter, L, feet. 

The difference between average base pressure and hub streamline 
base pressure divided by hub streamline tangential velocity 
head ; d imens lonles s . 


6 


"ml 

FF 


PL 


- 15 - 



RCUO 

DRCUO 


(P^/P,) 

T^O 

ROT 


IPX 

C. 


Angular momentum, RVg , on tip streamline; ft^/sec. 

Angular momentum difference between hub and tip, (RV - 
ft /sec. The distribution of RV is assumed 
linear with stream function between hub and tip. 

Nozzle exit pressure ratio at tip; dimensionless. 

Inlet temperature - upstream of rotor for rotor case; deg. R, 

0 if swirl is generated by a stator; 1 if swirl is generated 
by a rotor. 

1 if downstream flowfield printout is desired; 0 if not. 

Ratio of downstream hub streamline radius to nozzle exit tip 
radius; dimensionless. (See Figure 91 for values recommended 
from data.) 

Empirical constant determining mean streamline curvature 
gradient at nozzle exit station. (See Figure 92.) 


The following four inputs are used only for the rotor option but 
must have dummy values for stator option also; 


0 ^ Rotor tip solidity; dimensionless. 

a„ Rotor hub solidity, dimensionless, 

n 

Average rotor inlet absolute velocity; ft/sec, 

n Average rotor efficiency* 

R 

Output Definitions 


Output from the program is in three blocks. The first, which is 
optional (depending on IPX), describes the flowfield downstream of the 
nozzle exit. At each axial station, it lists in this order: 

V Hr- f 


The second output block is generated only for the rotor case and 
approximately describes the rotor blade environment, the rotor blade 
static pressure rise coefficient and "D" factor at each streamline. 
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At each streamline (starting at the hub), it lists in this order across 
the page: 

S»L« No> , R, ^ 

The third output block lists most of the input on the first line 
for reference. On the second line, it lists the overall performance 
parameters and other items of interest in this order: 

Cp, Cy, Crp, CpL, ^P2^ 

Where Z is the value at the downstream Station 2 and L is the final value 
of the length parameter. 
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EXPERIMENTAL INVESTIGATIONS 


TEST HARDWARE AND INSTRUMENTATION 


Base Pressure Models 

The base pressure model consisted of a 4 inch (10.16 centimeter) inner 
diameter stainless steel casing in which a 2 inch (5.08 centimeter) aluminum 
center body was positioned by one support strut as shown in Figure 5. To 
the inlet of the models a small bellmouth was attached. Vanes to generate 
axial or swirling flow shown in Figure 6 were fitted behind the center body 
with the swirl vanes designed to produce a comparable level of swirl as 
observed in the full scale testing of the LF336E (Reference 1). Inter- 
changeable inner and outer aluminum flowpaths could be positioned at the 
model exit plane. A total of 11 sets of different inner and outer flow- 
paths was designed with varying radius ratios and flow angles with a 
positive flow angle designated as radially outward. Table I contains perti- 
nent data on the exit flowpath and Table II gives the contour geometry of 
each set. Figure 7 shows a schematic of the test model cross section. 

For simplification, the model configurations will be hereafter desig- 
nated by the model number as given in Table I and a prefix of either an "A" 
or "S" will be used to identify axial or swirling flow, respectively. Thus, 
S5 is the test configuration designation for Model 5 tested with swirling 
flow and A5 is for Model 5 with axial flow. 

Test Facility 

The base pressure tests were performed at the probe calibration 
facility at the General Electric Company, Evendale, Ohio. This test 
facility has the capability of adapting the probe stand for small scale 
model testing. With the base pressure model positioned in the test stand 
as shown in Figure 8, total pressure ratios of slightly greater than 1.4 

and flow rates between 2-3 Ibs/sec (0.9-1. 4 Kg/sec) could be attained using 
the shop air system. 

The test facility contained instrumentation systems to provide constant 
temperature inlet air operation for each test and provided the actuation 
systems for the exit plane traversing probe in immersion, yaw and pitch. 
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Instrumentation 


Instrumentation on the base pressure model was kept at a minimum level 
with the major portion of the performance data being recorded by three port 
exhaust traversing probes. Eight static pressure taps were located 90^ from 
the center body support strut and equally spaced from the inlet of which 
the first five were used to record wall statics for upstream flow measure* 
ments. 

The three port exhaust traversing probe was mounted in a facility 
actuation system which had the capability of immersing, pitching and yawing 
the probe. For the base pressure tests, the probe was installed such that 
the probe static pressure ports could be nulled, i,e,,. probe aligned in the 
direction of the flow which yielded the magnitude and direction of the 
swirl angle. The data pressure null was indicated by a pressure balancing 
transducer with a null meter readout with each static pressure connected 
to each side. 

The exhaust traversing probe was mounted aft of the model such that 
the axial distance could be varied to at least four model diameters down- 
stream. 
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TESTING AND PROCEDURES 


Test Scope 

The main purpose of this test was to investigate the effects of various 
flow field conditions of annular jets on the base pressures of blunt. bases 
so that experimental data could be correlated with an analytical model in 
order to develop a useful tool in the future design of turbotip lift fan 
systems. The new analytical model, including empirical constants from the 
tests, could then be used to predict the base pressure levels of new lift 
fan designs, thus reflecting permissible or nonpermissible levels of base 
pressure or exit geometry. 

A total of 50 hours of testing was conducted consisting of three to 
four traverse probe axial positions with axial and swirling flows, variable 
pressure ratios, blockage at the exit plane to simulate louvers, and skewed 
swirl distributions. Table III gives a complete summary of the different 
configurations that were tested. 

Test Procedures 

Since a minimum of instrumentation was used for the tests and due to 
an unavailability of an adequate pressure recording system, the total and 
static pressures from the traversing probe and the wall statics were 
recorded manually from mercury manometers. Traverse data was recorded for 
approximately 20 radial locations from the tip to center of the base. 

Digital counters on the traverse probe yaw and immersion actuators were 
used for determination of the swirl angle and the probe radial location. 

Probe Calibration 

One three-port cobra probe was used to obtain a survey of the total 
and static pressures and the yaw (swirl) angle across the model exhaust 

plane. Figure 9 shows a schematic of the probe head with several pertinent 
dimensions. 

The probe was calibrated in the probe calibration stand over a range 
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of pitch angles (4*10** to -20**) and Mach numbers (0-6 and 0-8). The 
calibration was performed with the yaw angle kept constant (0°) since the 
actual model tests would be performed with the static pressures nulled^ 
i.e.^ the probe aligned in the direction of the flow. 


The recorded pressures were reduced into a coefficient form such 

that 


K1 = 


P - P 
TO a 

P - P 
TP SP 


where, 



- plenum total pressure 

- barometric pressure 

- total pressure indicated by probe 

- static pressure indicated by probe 


Also, a correction factor (K2) was included due to the effect of the flow 
angularity (pitch) , where 


K2 = 

Pgp (3 = 0“) 

The results of the probe calibration are shown in Figure 10. These 
calibration characteristics provided the basis for converting the measured 
pressure data into the appropriate corrected pressures, where 

Pj = P^p - (Ppp - Psp)(Kl/K2) 

Data Processing 

Data reduction of the recorded parameters was done by a time-sharing 
computer program which calculated and integrated various flow field param- 
eters. Further details of the procedures used for attaining overall model 
performance are given in Appendix II. 
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TEST RESULTS 


The results of the Base Pressure Tests are presented in the follow- 
ing section of this report. The discussion will be separated into major 
categories covering; 

• Model performance with axial flow 

• Model performance with swirling flow 

• Model performance with variable pressure ratio 

• Model performance with simulated louvers (blockage) 

• Model performance with skewed swirl distribution 

Performance With Axial Flow 

Data was recorded for the axial configurations at several traverse axial 
distances. Figures Uthrough 16 show the total and static pressure coefficient 
distribution for the closeup traverse position and Figures 17 and 22 show the 
same parameters but with the traverse plane being located 0.656 model 
diameters downstream. As indicated from these figures, the configurations 
with a positive flow angle, radially outward, exhibit a higher static pressure 
in the flow stream and a lower static pressure in the base region as those 
configurations with axial or less than axial flow angularity. 

Integration of the exhaust plane at several axial distances made it 
possible to locate the flow streamlines in the exhaust region of the models. 
Figures 23 through 28 show a graphical representation of the exhaust flow 
fields for all the models tested with axial flow. The 1007, streamline loca- 
tion was determined at the point where the integrated flow equalled the up- 
stream flow with the integration starting at the centerline and proceeding 
radially outward. The zero percent streamline was determined at the location 
at which the total pressure was equal to zero. The mixing region, as indi- 
cated on the figures, is the difference between the total integrated flow 
and the upstream flow. In this region, the moving jet interacts with the 
stationary surrounding air, transfers momentum to the stationary air, thus 
entraining flow. Figures 29 through 31 show the amount of flow entrainment 
for all the models with axial flow. 
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Figures 32 through 34 contain the wake' mixing loss characteristics for 
the models with axial flow with the mixing loss coefficient defined as, 


where 


- ^ ^TO ~ ~ ~ 

^TO ■ 


TO 

^Sm 

It 

O) 

o 


- upstream total pressure 

flow section static pressure 

- mass averaged total pressure 


- loss coefficient due to components of model (bellmouth, 
vanes, etc.) 

P - ambient pressure 

Figures 35 through 38 show the effects of model geometry on the average 
hub base pressure, the main flow thrust coefficient, the overall thrust co- 
efficient and on the flow coefficient for axial flow. The calculation pro- 
cedures for these parameters are discussed in Appendix II. 


Performance With Swirling Flow 


Figures 39 through 44 show the total and static pressure coefficients 
and swirl angle distributions for the set of models tested with swirling 
flow with the exhaust plane survey being made at the closeup traverse position. 
The same parameters for the same test conditions, but with the exhaust 
traverse plane being 0.656 model diameters downstream are presented in 
Figures 45 through 50. An immediate observation is the difference in the 
levels of the static pressure coefficients between configurations with axial 
and swirling flows. For the swirl cases, the static pressure coefficients 
are at least twice as low than those with axial flow. 


The graphical representations of the flow fields for the models with 
swirling flow are shown in Figures 51 through 56. The same calculation 
schemes were used to establish the 100 and the zero percent streamlines and 
the mixing regions as with the axial cases. Figures 57 through 59 contain 
the flow entrainment characteristics for the configurations tested with 
swirling flow. Again, an immediate observation shows that the flow entrain- 
ment characteristics of the configurations with swirl are 20 to 30 percent 
higher than the same configurations with axial flow. 
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The wake mixing loss characteristics for the swirl cases are shown in 
Figures 60 through 62 • The mixing loss coefficients are calculated in the 
same manner as those for the axial cases. 

The effects of model geometry on the level of the average base pressure 
and on the base pressure at the hub outer radius are shown in Figures 63 
through 65. As can be seen, swirling flow has a tremendous effect upon the 
levels of base pressures that were observed. Compared to the axial cases^ 
the swirl configurations are, on an average, four to six times lower, i.e., 
more negative. 

The effects of geometry on the main flow thrust and on the overall 
thrust coefficients are shown in Figures 66 and 67. Figure 68 contains 
the variation of the flow coefficient with model geometry. 

The unstable regions as indicated in Figures 63 through 68 were 
observed while testing the 0.4 radius ratio models at -10 and -20 degree 
flow angles. Intermittent, unstable conditions existed several times with 
these two configurations which were signified by changes in the audible 
noise frequency and a fluctuation of the total and static pressures in the 
exhaust flow. Inserting a tool or a hand into the exhaust stream removed 
the unstability. 

Performance With Variable Pressure Ratio 

The three models with zero degree flow angles were tested at two other 
pressure ratios, 1,2 and 1.4, with swirling flow to investigate the effects 
of variable pressure ratio upon the base pressures. Figures 69 through 71 
show the effects of pressure ratio on the total and static pressure coeffi- 
cients and on the swirl distributions. As indicated, the swirl angle and 
total pressure coefficients are unchanged by pressure ratio while the static 
pressure coefficients do show a slight change. Figures 72 and 73 show the 
effects of pressure ratio on the base pressure coefficients averaged over 
the base and at the base outer radius. The trends that are shown indicate 
that the average base pressure increases slightly with pressure ratio while 
the base pressure at the base outer radius exhibits the opposite trend. 
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The effects of pressure ratio upon the main flow and the overall thrust 
coefficients are presented in Figures 74 and 75* As shown, the thrust coeffi- 
cients increase with increasing pressure ratio. Figure 76 shows that 
increasing pressure ratio has no significant effect upon the measured flow 
coefficient. 

Performance With Blockage 

The possibility that the base pressures could be increased by the 
addition of a blockage system, simulating exhaust louvers, was investigated. 
The blockage system consisted of a wooden support frame containing eight 
rows of 0.0645 inch (0.1638 cm) outside diameter, stainless steel tubing. 

The tubing spacing was selected at 0.5 inch (1,27 cm) such that a 14 per- 
cent blockage existed. The blockage system was mounted 0.5 inch (1*27 cm) 
aft of the nozzle exit plane. Figure 77 shows photographs of the test setup 
with the blockage system mounted in position. 

The models with zero flow angles were tested with the blockage system 
Installed with and without swirling flow. Figures 78 through 80 show the 
effects of blockage with axial flow. The only result that the blockage 
produced was a small increase in the flow stream static pressure of approx- 
imately 5 percent for all three configurations. The influence of the 
blockage upon the base pressure and the flow coefficient is insignificant 
as shown in Figure 81^ 

Figures 82 through 84 contain the effects of blockage on the total and 
static pressure coefficients and on the swirl angle for those models tested 
with swirling flow. As with the axial flow configurations, a similar trend 
is observed in regard to the stream static pressure, that is, a slight 
increase in that static pressure. Figures 85 and 86 show the influence 
of the blockage and the base pressures and on the flow coefficinet. The 
blockage increased the base pressures approximately 10 percent and de- 
creased the flow coefficient by 2 percent. 
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Performance With Skewed Swirl Distribution 


Since the axial flow configurations exhibited a much lower level of 
base pressure as compared to the swirling flow cases, the possibility that 
a skewed swirl distribution across the flow stream would improve the base 
pressures was investigated. 

The axial vane assembly used to produce axial flow was reworked such 
that a similar level of swirl as observed with the swirl vane assembly 
existed at the tip of the vanes. At the hub, no rework was done so that 
axial flow would exist there, thus intending to produce comparable levels 
of base pressure as observed in the axial flow cases. 

Figures 87 through 89 show the total and static pressure coefficients 
and swirl distributions for the Models 2, 6 and 10 with zero flow angle. 

As shown, the levels of base pressure are between those observed for the 
nxial and swirl cases. Table IV contains a performance comparison for 
Models 2, 6 and 10 with the three types of flow fields tested. These 
results show that for the skewed swirl distributions, base pressure, 
thrust and flow coefficient all increased from the observed level with 
complete swirling flow, the largest increase occurring in the thrust and 
flow coefficients. 
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APPLICATION OF ANALYTICAL METHODS 


The analytical model developed for representation of the flow around 
blunt bases requires an evaluation of numerous empirical constants for 
completion of the mathematical model. The four most significant constants 
are : 

FF - Effective friction or loss factor representative of losses 

associated with the boundary between the main and wake flows. 

C 2 - The area occupied by the wake flow field at axial locations 

far downstream of the hub base. This constant represents the 
ratio of the wake flow diameter to the nozzle tip diameter. 

Cp^- The coefficient representative of the pressure gradients in the 

wake flow at the plane of the hub base. The parameter establishes 
the integrated average base pressure in terms of the pressures at 
the outer radius of the base. 

Cg - A coefficient that controls the initial rate of change of curva- 
ture of the wake flow. 

The tests of the base pressure models with swirling flow provided 
the data necessary for evaluation of these constants. The initial step 
in the analysis was to determine the appropriate flow field parameters in 
a format compatible with the analytical model. The inlet flow to the 
model is represented as an average total pressure ratio, an air total 
temperature and the angular momentum distribution. The total pressure and 
temperature inputs were determined through test measurements and were 
assumed constant for all test models. The angular momentum in the analysis 
is represented by a tip momentum with a linear variation between the hub 
and tip. The test measurements were used to evaluate the required momentum 
distributions by fitting the measured profiles. The resulting variation of 
tip momentum and hub-to-tip momentum difference is given in Figure 90 for 
the range to test radius ratio. The total pressure ratio and temperature 
used during the analysis were 1.263 and 535® R (297® K) , respectively. 
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Only one of the four empirical constants can be obtained directly from 
the test measurements. The wake area coefficient, C 2 , can be obtained by 
referring to the test flow streamline distributions. Using the most down- 
stream traverse location, the radius of the wake region was obtained and is 
summarized in Figure 91. The wake size appears to be Independent of both 
model radius ratio and radial flow angle. An average value of 0.48 was 
selected for the coefficient, C 2 ^ for evaluation of the test data using the 
mathematical model. A constant value of the wake size appears reasonable 
for the test model since the model employed one set of turning vanes up- 
stream of the contoured nozzle and exhaust plane. The swirl vane had a 
hub- to- tip ratio of 0.5 with a tip radius of 4 inches (10.16 cm). Through 
observations of the downstream flow profiles, it appears that all models 
exhibit similar swirl and pressure distributions. Apparently, the flow 
field, as established at the plane of the swirling vanes, determines the 
downstream flow conditions and is not affected by the internal flowpath 
of the actual test model. This consistent flow field pattern can be 
expected since both axial and tangential momentum are conserved except for 
friction and mixing losses. 

With one of the constants established, the problem was to determine 
the three remaining constants using the model test data. Through a process 
of trial and error, the friction factor was established at a level of 0.04 
and the values of Cg and Cp^ were determined to be a function of radius 
ratio and radial flow angle as shown in Figure 92. The justification for 
these characteristics is based on agreement of predicted and measured hub 
base pressure levels, A comparison of the hub base pressure coefficients 
as obtained from tests and theoretical analysis is given in Figures 93 
through 95. 

The agreement of test and theory is fair considering the small 
sampling of test configurations and the complexity of the analysis. A 
comparison of test and analytically derived nozzle velocity and flow co- 
efficients shows similar agreement. The agreement of velocity coefficients 
is of prime importance since it includes the combined effects of swirl 
^^®t^^^utions in addition to stream and base pressure levels. The agreement 
between test and theory is within 2 percent. 
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The empirical constants^ as derived for comparison of configurations 
with swirling flow, were then used to evaluate hub base pressure levels 
during axial flow* Without swirling flow, the value of the constant C 2 
should be very small since the jet wake closes on Itself at the hub base 
centerline* Likewise, the base pressure coefficient, would be very 

small since the stream static pressures are uniform throughout the separated 
base regions. Using these assumptions, the comparison of test and analysis 
as shown in Figures 96 through 98 is possible. Again, the comparison shows 
fair agreement with the largest differences occurring for the 0.4 radius 
ratio configuration* 

The configurations tested with the skewed swirl distribution and 
model E were compared x^Tith the analytical model using the empirical con- 
stants that were determined previously. Table V contains the comparison 
of test and theory* As shown, the theory predicts the test results 
reasonably well. This final comparison verifies the selection of the 
empirical constants derived from the 20 configurations previously tested 
and also confirms the validity of the analytical model. 

Based on the comparison of test and analytical parameters such as 
velocity, flow and pressure coefficients, the analytical model appears to 
be an adequate representation of the exhaust flow downstream of annular 
nozzles with blunt bases. Considerable additional testing, including 
numerous models with variable swirl levels, along with further refinement 
of the analytical model, would be required to improve the levels of data 
correlation. 


- 29 - 



CONCLUSIONS 


An analytical model was developed to predict the base pressures 
established by axial or swirling flow through annular jets having 
blunt bases. Configurations consisting of various flow angularities, 
radius ratios and swirl distribution can be modeled by the computer 
program generated from the developed theory. Over 50 hours of 
scale model testing were conducted to assist in the selection of 
several empirical constants required by the analytical model. The 
significant conclusions which were derived from the results of this 
investigation are as follows: 

• Three separate factors have been recognized as major contributors 
to the pressure at the base of an annular nozzle, i.e., the shear 
forces between the high velocity main stream flow and the rela- 
tively stagnant air mass composing the centerbody wake, the 
meridonal curvature of the flow streamlines in the vicinity of the 
nozzle exit, and the radial pressure gradient due to the swirl 
velocity component. 

• Major differences in the flowfields between axial flow and 
swirling flow were observed. With axial flow, base pressure 
coefficients of -0.2 to 0,25 were measured while with swirling 
flow, base pressure coefficients of -0.4 to -0.8 were not 
uncommon. 

• Base pressure levels are unaffected by changes in pressure ratio 
and only slightly affected by blockage downstream of the base. 

• The agreement of test and theory is fair considering the small 
sampling of test configurations and the complexity of the analysis. 
The analytical model is an adequate representation of the exhaust 
flow downstream of annular nozzles with blunt bases based on the 
comparison of test and analytical parameters. 
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• Additional testing is required to refine the analytical model 
in order to improve the correlation between the theory and the 
test. 
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NOMENCLATURE 


Symbol 

Definition 

a 

n 

Acceleration normal to a streamline 

A 

Constant, 1 - R * 
m2 

A 

m 

Model measuring section effective 
flow area 

A 

n 

Nozzle exit annular area 

AP 

Surface area of wake mixing zone 

& 

Constant, A + L tan 6 , 

ml 

c 

m 

Mean streamline curvature 

c 

Constant, (-A + 3B + C^)/6 

S 

Nozzle flow coefficient, M/M^ 

S 

Specific heat at constant pressure 

Sl 

Empirical coefficient 
(\ - Pgi)/tl/2 »i(RV,)\j 


Rotor static pressure rise coefficient 

<^2 - 

CPS 

Local static pressure coefficient, 
‘■■s - 

CPT 

Local total pressure coefficient, 

O-T - 5^)/(P„ - P.) 

Si 

Hub base pressure coefficient at outer 
radius, (P^ - - P^) 

S 2 

Average hub base pressure coefficient, 

- V/<^i - 

s 

Nozzle thrust coefficient, F/F^ 

s 

Nozzle velocity coefficient, C_/C„ 

T r 
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Untts 

ft/sec^ (m/sec^) 
ft (m) 

ft^ (jnh 

ft^ Cm^) 
ft^ (m^) 
ft (m) 
l/ft (1/m) 
ft (tn) 

BTU/lb-"R 

(joule/kg-*K) 



Definition 


Units 


Symbol 

Co 


d 

DF 

DRCUO 


®1 


P 


FF 


g 

IPT 

J 

K1 

K2 

L 

M 


Ratio of downstream to nozzle tip 
radius, Rt2^®T1 

Empirical constant 

Total differential 

Rotor diffusion factor 

Angular momentum difference between hub 
and tip, CRVq)jj - CRVq)^ 

Iteration error function in radius ratio 

Iteration error function In axial force 
balance 

Total thrust 

Ideal thrust 

Wake shear stress coefficient, 

F /(1/2 p V 2) 
b m 

Gravitational constant 

Print control parameter, 1 if flow*- 
field print is desired, 0 if not 

Work equivalent of heat 

Probe calibration coefficient, 

(^T “ ^s^^^**TP " ^SP^ 

Probe calibration factor. 

Length parameter 
Mass flow rate 

Ideal mass flow rate 


n 


Distance normal to streamline 


ft Cm) 

2 2 
ft / sec C® /sec) 


lb (kN) 
lb (kN) 

2 2 
ft/sec (m/sec ) 

ft-lb/BTU 
(N-m/ joule) 


ft (m) 

lb-sec/ ft 
(kg-sec/m) 

Ib-sec/ft 

(kg-sec/m) 

ft (m) 


- 33 - 



Symbol 

Definition 

Units 

P 

Static pressure 

Ib/ft^ (kN/m^) 

P 

a 

Ambient pressure 

Ib/ft^ CkN/m^) 


Area averaged hub base pressure 

Ib/ft^ (fcN/m^) 


Corrected probe static pressure 

Ib/ft^ (VN/mh 

^Sm 

Model flow section static pressure 

Ib/ft^ (kN/m^) 

^SP 

Probe indicated static pressure 

Ib/ft^ (kN/m^) 


Total or stagnation pressure 

Ib/ft^ CkN/m^) 


Mean total pressure at nozzle exit 
based on ideal velocity, 

Ib/ft^ (kN/m^) 

^TP 

Probe indicated total pressure 

Ib/ft^ (kN/m^) 

^TR 

Relative total pressure 

Ib/ft^ (kN/m^) 

Q 

Dimensionless indicator (0, stator 
case; 1, rotor case) 


r 

Streamline radius of curvature 

ft (m) 

r 

m 

Mean streamline radius of curvature 

ft (m) 

RCUO 

Angular momentum at tip streamline. 

ft^/sec (m^/sec) 

ROT 

Control parameter, 0 If swirl is 
generated in stator, 1 If not. 


R 

Radius 

ft (m) 


Hub radius 

ft (m) 

«m 

Mean streamline radius 

ft (m) 

«T 

Tip radius 

ft (m) 

S 

Distance along a streamline 

ft (m) 

T 

Static temperature 

®R (”K) 

T„, 

T 

Total temperature 

'’R (*“K) 
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Symbol 

Definition 

Units 


U 

Wheel speed 

ft /sec 

(m/sec) 

VIN 

Average rotor Inlet absolute velocity 

ft/ sec 

(m/ sec) 

\ 

Hass averaged ideal jet velocity 

ft/sec 

(m/sec) 


Meridinal velocity component 

ft /sec 

(m/ sec) 

\ 

Total velocity component at tip 

ft/sec 

(m/sec) 


Axial velocity component 

ft/ sec 

(m/sec) 

■ 

Tangential velocity component 

ft/sec 

(m/sec) 

w 

Integrated airflow from wake 
centerline 

lb /sec 

(kg/ sec) 


Rotor inlet relative velocity 

ft/sec 

(m/ sec) 

W2 

Rotor exit relative velocity 

ft/sec 

(m/ sec) 

T 

Total integrated airflow at nozzle 
exit plane 

lb/ sec 

(kg/ sec) 

XIMM 

Probe immersion 

ft (m) 


Z 

Axial distance downstream of 
nozzle plane 

ft Cm) 


a 

Yaw (swirl) angle 

deg (deg) 

P 

Streamline merldonal flow angle 

deg (deg) 


Mean streamline merldonal flow 
angle 

deg (deg) 

Y 

Ratio of specific heats 



3 

Partial differential 




Average rotor efficiency 



0 

Tangential flow angle 

deg (deg) 

X 

Angular momentum distribution factor, 
-I(RVq)^ - (RVg)jj]/ii 
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Symbol 

P 


Definition 


Mass density 


Ideal mass density 


a 

01 


Mass density based on stagnation 
conditions 

Solidity 

Stream function 

Wake mixing loss coefficient, 

CP„ - V/^n - 

Nozzle ^del Internal loss coefficient. 


Units 

Ib-sec^/ft^ 

(kg-sec^/m^) 

Ib-sec^/ft^ 

(kg-sec^/m^) 

Ib-sec^/ft^ 

(kg-sec'^/m^) 
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Subscripts 

0 Upstream or reference station 


1 Nozzle exit station, or rotor Inlet station 

2’ Downstream station, or rotor exit station 

H Hub of flow field 

m Model flow or measuring station 

T Tip of flow field 
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APPENDIX I 


DERIVATION OF EQUATIONS OF EQUILIBRIUM NORMAL TO A STREAMLINE 
Angular Momentum Distribution 

It is assumed in this analysis that the angular momentum distri- 
bution is known at the nozzle exit station and that it can be adequately 
represented as a linear function of the stream function Y as follows: 

RVq = (RV^)^ + X M (1 - Y) (ID 

where 

, - I(RVj)t - (RVg)g]/ft (12) 

Assuming that angular momentum is conserved along any streamline, 
the above two equations can be applied at any point downstream of the 
nozzle exit station as well, 

Pressure-Temperature-Density Relations 

Assuming uniform stagnation pressure and temperature at the upstream 
Station 0 and an isentropic process between Station 0 and the nozzle 
exit Station 1, we write 

P = P fT /T (y~1) 

which will apply to all streamlines at Station 1. Further, if we know 
the stagnation pressure loss between Station 1 and any downstream 
Station^ we may write 

P/^TO = (13) 

which will apply to any point downstream of the nozzle exit. Similarly, 
we may write for the density: 

p/Pto “ 

Notice that the above three equations apply even though there may 
be energy addition (as by a fan rotor) between Stations 0 and 1 as long 
as the process is isentropic. 
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Energy Relations 
Case 1: 

If no energy is added between Stations 0 and 1, then conservation of 
energy along a streamline results in 

’"/■'to ■ 1 J s V 

Case 2; 

If the angular momentum at the nozzle exit is entirely the result of 
isentropic energy addition by a fan rotor^ then conservation of energy 
along a streamline results in 

T/Tto - 1 + [Vjj2(RVj/(RV,)j) - v/ - V^2]/(2g J Cj Tjp) 
where the energy added at the rotor tip is 
AT^ = V^^/(2g JCp) 

We may combine the above two cases into one equation 

T/T^o ■ 1 J S Tto> 

by defining 

Q = 0 for no energy addition (stator case) 

Q = 1 for rotor energy addition (rotor case) 

Since Vg = RVg/R the above equation may be written as 

, (RV„)2 

T/Tq “ ^ 1 \ y(2g J Cp T^g) (15) 

R 

Continuity 

By noting that 

dR = -dn cos fl 

m 

we may express the continuity equation along any stream tube as follows: 

dY = 2 tt R p V dR/M cos S 
or 

SR _ M cos S 
5Y ^ — !2- 

2t 7 p R (16) 
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Momentum Normal to a Streamline 


The acceleration normal to a streamline of an element of fluid is 
given by 

The mass of this element is given by 


dm = P R de dS dn 

The normal pressure force which balances the momentum force is given by 
dF = R de dS dP 


Then since 


dF = a dm 
n 


we may combine the above three equations to write 

g - P (V„'/r - Vg' CO, VW 

The derivative of stream function may be expressed as 

dV « • 

jhT “ -2^7 R p V /M 
on ^ m 


and the last two equations may be continued to yield 

II = (Vg^ COS - V^Vr)/(2tt R V^) (17) 

This equation gives the pressure gradient normal to a streamline required 
to balance the momentum forces normal to the streamline. 

Combined Equation 

Equations II through 17 can be combined to eliminate pressure as a 
variable in Equation 17 and yield a new form for the momentum equation. 

This will then eliminate the need for Equation 13 in subsequent calculations. 
First, we differentiate Equation 13 with respect to the stream function to 
get 



where we have assumed that the total pressure loss factor is constant in 
a direction normal to the streamlines so that = 0. Using the 

perfect gas equation of state and Equation I4 in the above result gives; 


SY 


9(t/t ) 

g J Cp p — 


( 18 ) 


Next^ we differentiate Equation 15 with respect to stream function and use 
Equation II in the result to get 


MT/Tg) . ^ ; X Q _ 2 r (rv,) ifK - (RVg)^ || . ^ , lil 

aY (R Vq)^ 

(2g J Cp T^q) (19) 


Using Equation 19 in Equation 18 to eliminate the temperature derivative 
and then replacing the left-hand side of Equation 17 with the result gives 


P C ^T1 


M XQ 


(RVe)^ 


2r (RVq) M X - (RVq) 


2 M 
aY 


. 2 V 1 = 

m QY ^ 


vj- cos 3 /R - V ^/r 
y mm 


TT R V, 


m 


Use of Equation 16 in the above and then solving for 9 finally gives 

a new form of the momentum equation as follows: 


S V 


m 


9Y 


M 

2n p R r 


+ 


t 0 


(RVj) ] 


R 


( 110 ) 


This equation gives the meridional velocity gradient required to balance 
the momentum forces normal to a streamline • 
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APPENDIX II 


DATA ANALYSIS PROCEDURES 


The exhaust flow of the various models was surveyed by radial 
traverses using the 3 element probe^ The flow surveys for the various 
configurations consisted of approximately 20 discrete data points 
approximately 0.1 inch (0.254 centimeters) apart. The outer radial 
traverse boundary was established at the location where the total 
pressure was just equal to zero. The inner radial boundary was 
selected on the centerline of the particular model being surveyed. 

The following discussion presents the calculation procedures used 
for analysis of the exhaust probe data. The calculations were performed 
using a time-sharing computer program. 


The following data obtained from the exhaust probe was used as 
input to the data reduction program: 

• XIMM, Probe Immersion 

# Pj^ Total Pressure (corrected for probe calibration) 

• Pg, Static Pressure (corrected for probe calibration) 

# Yaw (Swirl) Angle 


Additional fixed data input based on average conditions during 
the complete traverse cycle where: 

• Total or Plenum Temperature 

• P^, Barometric Pressure 

a 

• Model Flow Section Static Pressure 

• 6, Model Flow Angle 

• Upstream or Plenum Total Pressure 


Several overall performance parameters were initially calculated 
as follows: 

p 

1/2 


• Flow Section Mach Number, M = [( — ^ 

Y-1 


- !)]• 


• Flow Section Static Temperature, = T„/(l = M 

Sm T V ttt •' 


Flow Section Velocity, V - M (v g R ) 

m m Sm 


T ' Y 
1/2 
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• Flow Section Density, p = P_ /g R T„ 

* m Sni Sin. 

• Mass Flow, M = p V A 

m m m 

• Nozzle Exit Total Pressure, = P^^ -- (P - P ) 

> T1 TO 0 ^ TO Sm^ 

p 

• Ideal Mach Number, = [( ^ 

3 - 

• static Temperature, Tg = T^q/(1 + 

• Ideal Velocity, V, = M, (y g P 

IX o 

V 

• Ideal Mass Flow, M. = (P g R T„)V, 

X a ® S i 

9 

• Ideal Thrust, F^ = 


At each discrete point or immersion, the following parameters 
were calculated: 

Y-1 

P Y 

t Absolute Mach Number, M = [ ( — - 1)]^^^ 

Y-1 Pg 

• Total Pressure Coefficient, CPT = (P - P )/(Y - P ) 

0 Static Pressure Coefficient, CPS = (P - P )/(P - P ) 

S a XI a 

• Static Temperature, T = T^/ (1 + M^) 

• Absolute Velocity, V = M(y g R T)^^^ 


• Flow Per Unit Area, W/A = (Pg/RT) V cos a cos 6 

• Thrust Per Unit Area, F/A = W/A V cos a cos S + (P - P ) 

S a 

The total model exhaust flow was integrated to obtain the follow- 
ing performance; 


• Airflow, W = Z (W/A) AA 

• Total Thrust, F = E (F/A) AA 

• Mass Average Total Pressure, P^^ = (zp^ (W/A) AA)/W 

• Hub Base Pressure, P^^ = J:(Pgg) AA/A 

The following overall performance parameters were then calculated; 

• • 

• Flow Coefficient, C„ = M/M, 

F i 

• Thrust Coefficient, C^ = F/F^ 


- 44 - 



• Base Pressure Coefficient at Hub Exit Radius, C. 


PI 


P - P 
HI a 

- P 

Tl a 


• Average Base Pressure Coefficient Across Base, C 


P2 


Pu - P 

H a 

Pn,i - P 
Tl a 


The flow streamline locations were determined as follows: During 

the flow integration, the radius and the resulting integrated flow were 
stored in an array throughout the entire flow stream; then, for 10 to 
100 percent flow, the appropriate radius was determined from the stored 
data by linear interpolation. 
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APPENDIX III 


COMPUTER PROGRAM LISTING 


1*#RUNH * =HBASE(»OGO,CORE=30) 

100 COMMOH PRATIO,PAM,ROTO,PLOSS,RKOl,RCIJ.PI,RHOG,V0,P0,MF,RCU0,Un 

I10i ,GI,G2.C3,CA.PT1NT,OMF.T0,ROI.SQU,SCLH,IND,VIN,EPF 
120 REAL LAM,MF,DMF(m,MFI,L,RlP<A),R0P(4),ZlP(A>,Z0P(4),BKP<4) ,ZMP(4) 

130 INTEGER ROT 

140 DIMENSION RPPU) ,X<3) ,Y(5) ,UI ( 12) ,VI2< 12) ,E1 (5> ,E2(3) ,RCU(12) , 

1504 PRATI0C12),R<18) 

160 1 PRINT i00 } PRINT 100 

170 PRINT, R/R,BC,KF,L,CPL,RCU0,DRCU0 | READ,RB,B0,FF,L,CPL,RCU0,DRCU0 

180 IF<RR*EQ.0.0) GO TO 999 

190 PRINT,“P/P,TIH,ROT,IPT,C2,C9,C10,SOLT,SOLH,VIH,EFF" 

191 READ, V0, TIN, not, IPT,C2,C9,C 10, SOLI ,SOLH,VIN,EFF 

200 6=32. 17405; I NO=0 

210 GAN=1.4 { GUGAM/CGAM*!) ; G2sl/(GAM-1> 

211 T0=TIN»(l+<V0^<l/QI)-l)+(l/EFF-l)«ROT) 

212 RG=53.35 { RJ=778.i6 } C3=2X:Q4‘Gl«nG4<T0 

215 V0=SCRT(C3*<V0*U/Gl)-l)=mN/T0> 

220 ROTO=0 

230 1F(RCU0.N£.0.) ROTOsROT/RCU0 

240 C4=V0*2*.5*ROTO 

250 PAMsl4,696«t44 

260 P0SPAM 

270 IF<ROT.EQ.0> P0=PAM*<l-V0*2/C5)‘(«Gl> 

280 RHO0:P0/G/RG/T0 

290 PI=3. 14159 I B0=B£M>P1/I80 

310 1TR=I ; MAX1T:10 

320 C1 = P1*(1-RR‘'2)/(1-FRR‘^> 

322 VJ1AVs(V0*2-nCUa‘*2*<t+RR“2)/2)*.5 

324 TRATI 0= 1+ (V0"2« (RCU0+ .5*DRCU0)*ROTO-(RCU0+.5*DRCU0) ‘^-VWAV2)/C3 

326 RHO=RHO0*TRATIO'‘G2 

330 MF=2*Cl*RHO*VMAV*COS(B0) 

335 MF= ,9^MF 

340 RHl=RR*<2/<l+RR'‘2))*.5 

350 RH2=C2*(2/<1+RR*2))'‘.5 ; R«2sl 

360 V2:.99«V0 

370 DO 4 1=1,12 

380 RCU(I)=RCU0^DRCU0*(12*1}/10 

385 RCU(l)=RCU0fDRCU0 | RCU(2> sRCU&F.Be^ORCUO 

390 PRATIO(I) = <1+V0'‘2><‘RCU<I)*ROTO/C3>‘‘6UP0/PAM 

400 4 CONTINUE 

401 Z=2 
404* 

405 PLOSS=U*V2'^/<C3+V0“2*ROT))*(-Ql>/PRATIOU2) 

430 5 CALL ROTFLOU,RH2,0..0,,RT2,RH2,V2,FP2,FPH2.F«2) 

490 A=l-RH2 f . 

500 B=A4-L*TAN(B0) 

505 C=(-A+3*B+C9)/6 

510 1=1 ; DZMsZ/10 ; Z=0 

520 V2=V0 I AP=0 ; FFZ=0 * PFZ=0 ; PLOSS=l t 11=0 

525 BMAX=ABS(B0) 

530 10 £P=EXP(*Z/L> 

540 R«=<A+B*Z/L*C*Ziei2/L/U»EP+RM2 

550 R«P=**CA-BH<B-2*C>*Z/L+C*Z*Z/L/L)*EP/L 

560 RI1PP=<A-2*B+2*C+<B-4*C>*Z/L+C*Z*Z/L/L)*£P/L/L 

570 BM=ATAH(R«P) 

575 IF(ABS<3M>,GI,BHAX) B«AX=A8S(BH> 
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APPENDIX III Ccont'd) 


580 
590 
660 
602 & 
610 
620 
630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
750 
740 
750 
760 
770 
780 
790 
800 
810 
620 
830 
840 
850 
860 
870 
860 
890 
900 
910 
920 
930 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1045 
1050 
1060 
106U 
1065 
1070 
1080 
.1090 
1100 


C»=RMPF/( l+RMP-*2) *1,5 

CALL R0TFL0<5,Rt1,BF],C«,RT,RH,V2,VMT,VKK,PH> 

I F( <AB8( (RM2-RM)/(Rt52-l ) ) .GT. *01 ) .OR. (ABSCBM) .67. .2*B«AX> ) 

GO TO 13 

IFCIT.EQ.DGO TO 13 

Z^Z-DZFpKRflE+.BlstftBSd-RMai^ABSCRML'RKBl/CRnL-RKai-RMl/CRML-Rri) 

IT=1 

GO TO 10 

13 IFaPT-HNO.EQ.2) PRINT 2,Z,RMtB«*l80/PI,Cp},RT,RH,V2,VMT,VI4H,PH-PAM 
2 FORMATUH *F6,3,F7,5,F6.l ,F9.4,2F7.3,3F7.1,F8.0) 
ZHsZ+CRW-RH)*SIIKBH) 

ZTsZ-(RT-RM)*SIN<B«) 

IFU.FQ.l) GO TO 12 
BHsATAH((RH-RHL)/<ZH-ZHL)> 

0AP=PI * <RH+RHL )* <ZH-ZHL> /COS CBH) 

AP=AP+DAP 

DFF=.5^RH0li'FF*DAP*(VMH*2+V«HL*2)/2 

PL03S=PL0SSX><1-DFF/PTIHT) 

V2=<<C34-V0''2»>ROT>«<1-(PLOSS«PRATIO(12))*(-1/G1)>)*.5 

FFZ=FFZ+DFF*COS<BH) 

PFZ=PFZ+(<PH+PHL>/2“PA«)*DAP*SIN(BH) 

IF(aT.EQ.l>«OR.(I.6T,50))6O TO 20 
12 RHL=RH J RML=.RP1 J ZHLSZH > VKKLsVMH f PHL=PK 
IFd.GT.S) GO TO 15 

RIP<I)=RKL t ROP<l)sRT | ZlP(l>sZHL i ZOPCI)=ZT t ZWPUlsZ 
RPPCDsRM 
15 I = I-H t ZSZ+DZH 
GO TO 10 

20 ZIP(4)£0 j ZOP(4)s0 
CALL INT<RIP,Z1P) 

CALL 1KT(R0P,Z0P) 

RRT2RIP<4)/R0P<4) 

CALL IHT(Zi1P,ZIP) 

CALL IMKRPP.ZIP) 

RHURIP(4) t RMI=RPP(4) t ZW|sZMP<4) 

EPsEXP(-Zm/L> 

RMPs~(A-B+<D“2«C)*ZMl/L4C*ZNKZ«l/L/L)<^EP/L 

RMPP=<A-2*8+2*C-KB-4»s=C>*ZMl/L4-C*zm>l<zm/L/L)*EP/L/L 

BMl=ATAM<RMP) 

CH1=RMPP/(1+RMP*2)*1.5 

CALL RCTFLO(2,Rm.B«i,CMl.RTl,RHl.V0,FPl,FPHl,Fm> 

CALL ROTFLOU,RK2,0..0..RT2,RH2,W2,FF2,FPH2,FM2> 

ERRl=<RRT«RR)/RR 

FPHXiCPL*RHOl/2*PI«(RCU0+LAM*HF)*2 

FPHY=C10*FPHX 

ERR2=(FP2+FPH2-FFi2-FPl+FPHY-FPHI+FMl-FPHX)/FWl 
I F< ( (ABS( ERR 1 ) «LT. .003) . AND. (ABS (ERH2 ) .L7 . .003 } > 
«AND«((INO.EQ.U.OR.(1PT.EQ.D>>> GO TO 30 
1F((ABS(ERR1).LT..083).AND.(ABS(ERR2>.LT..003>> IND=1 
IFdTR.EQ.MAXiT) GO TO 26 
23 DO 40 1=1,2 
X(4-1)=X(3-I) 

Y(4-I)=Y(3-I) 


BEPRODUC1EE.ITY OP TE® 
PAGE IS POOR 
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APPENDIX III Ccont'd) 


40 


A2 


4t 

A2 

2.6 


30 


32 


34 


U10 
1120 
1130 
1140 
1130 
1135 
1160 
1165 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1350 
1360 
1370 
1380 
1350 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
14806 
1490 
1500 
1510 
15156 
1520 
1530 

1540 100 
1550 101 
1560 102 
1570 105 
1580 999 
1590* 
1600 
1610 
1620 
1630 


El(4-I)=El(3-l) 

E2(4“I>sE2(5-I) 

X(I)=KF ? Y(I)=L ; EJUlsERRi ; E2a) = ERR2 
IFaiR.EQ.l) GO TO 41 } lFCiTR.EEl,2) GO TO 42 
CALL IHT2<X,Y,El,E2,nF^L) 

PRIMI,X(l>,Y(l),ERRi,ERR2 
ITR=ITR+l I LAR:bRCU8/HF 
IFdKO.EQ. 1) PRINT 10^ 

GO TO 5 

L=U(l-*05*ABS<ERR2)/ERR2) ; GO TO 45 

MFrNF«U~RR" 2 )/<l+RK'‘ 2 )/(l“ftRT'' 2 >*<l+RRT* 2 ) I GO TO 45 

PRINT 101,HAXIT,ERR1,ERR2 

PRINT, ADDITIONAL ITERATIONS* j READ.NIT 

IF(NIT,Ee.G) GO TO 50 

MAXlIsMAXIT+NIT f GO TO 23 

NFI--0 s F1=0 

KHOI rR HO0* (PAM/P0) 1 /GAM) 

DO 32 Isl,l2 

CONTINUE^ 

C6s,2*VIU)-f »5*VI<2)+.5*VI(3> 

DO 34 l£3,ll 

C6=C6+<VI(I)-fVl(I+l))/2 

CONTINUE 

KFIspI*RHOI«<ROP(4>*ROP(4)-RIP<4)*RlP(4))*C6/10 

FIsMFl*C6/10 

CF=«F/MFl 

CTs (FM2-FP2-FPH2-FPHY)/FI 
CV=CT/CF 
VIE22<FI/«FI)*2 
C5 sC 3+V0“2*RCU(7)*ROTO 
DELPTs ( < 1*" VIB2/C5 )*<*Gl )- 1 )*PAM 
CPU-FPHl/0ELPT/PI/RHr2 
CP2s CP I* (FPH 1+FPHX )/FPH I 
PRINT 400 


PRINT, 
VI 


T0 


RR 


ROT 


60 

ITR’ 


FF 


CPL 


RCU0 DRCU0 


[a2,RR,B0*160/PI.FF,CPL,RCU0,ORCU0,C6/10,T0,ROI,lTR 


CF 


CV 


CT 


CPI 


CP2 


PRINT 
PRINT 100' 

PRINT,' 

MF FI’ 

S^TO 

). 

, "after ",I2,* iterations err is ",F6.3. 
,F7.5,F7.1,F8.4,F7.3,4F8.1,I5,I4) 
,5F7.3,2F7*2,F7.5,F7.2,F7.0) 


FOR mat UH 
FORMAKIK 
FGRMATUH 
FORMATUH 
STOP j END 


AP 


ERR2: 


SUBROUTINE lNTtR,Z) 

DIMENSION R(4>,Z(4> 

C2s (R (El'-R U ) ) /(z ( 2 ).z (I } ).c3* (Z (2>+Z U > ) 


L 


",F6.3) 


Z<2)) 
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APPENDIX III Ccont’d) 


1640 

1630 

1660 

1670* 

1680 

1690 

17006 

1710 

1712 

1715 

1720 

1750 

1735 

1740 

1742 

1744 

1746 

1750 

1760 

1770 

1780 

1VS0 

1800 

1810 

1820 

1C30 

1840 

1850 

1860 

1870 

18S0 

1890& 

1900 

1910 

1915 

1520 


CI=RU)»C2»2<l)"C54eZ<l)“2 
R<4) = C{+C2*Z<4)+C3#ZC4)'^2 
RETURN { END 


3 


5 


SUBROUTINE R0TFL0(KR,RH,Bi1,CM,RT,fiK,V2.PR 1 .PR2.PR5) 

COMMON PRATlO,PAri,ROIO,PLO2S,RHOl,SCU,PI,RHO0,Vf5,P0,MF,RCU8»LAM 
,Gl,G2*e5,C4,PnHT,DMP,T0,ROT,SOLT,SCtH,lND,VIN,EFF 
REAL LAM,VB(12),RU2)jp(12),RCUU2),PfiATI0(12),ttF,DMF<ll)»VX(l2) 
INTEGER ROT 

DATA UX/l.,l. Jr«l.,l..U,l*,l.,l.,l.J, I ,/ 

KSr7 I PLOS=PLOr>S 
IFCKR.EQ.DKSsi 


IF(:{R.EQ.2> PLO3=1.0 
K>0 $ PR1=0 : PR3;0 
DO 3 1=1,11 
DMFa> = .l#nF 

DHF<I>=.04*MF 8 DMF<2)=.06*MF 
VM(KS)=V2 t fi<l)=HH S R(7)=RH 
IF<AQS(CM).LT,2«E-18> CM=2.E-18 
DO 10 J=KS,11 

TRATlO=i+(V0'‘2=iRCU(I)4ROTO-(RCUa)/Ra>)'*2-V«a)‘‘2)/C3 

RHO=RHO0*IRATIO'‘G2*PLOS 

Fi s 1 /CM-l- (RM-R { I > > /C OS<BM ) 

F2-SIN<BM)*2 5 F3=RII)*2-RM*2 
F4=R(I)'‘24-<RM*SU1(BM))^2 

CVsI/F 14-F2*F5/F4''1,5 ; 1F(ABS<CV>«LT.2.E*18> CV=2.E-18 
F5=CV/2/PI/HHO/RtI) 

F6=LAM/Va<I )*(RCU<1 ) /R (1 )^2-C4 ) 

0VM1=F5^-F6 

DR 1 = . 5/PI /RHO/VMU ) /S (I )«C OS <BM > 

D|«A T 1 0= (-C4<=L AH+R CU Cl ) /R (I ) *2*( LAM4-RC U (I ) /R (I )*DR 1) -VM< I ) CD Vfl I ) 
DRH0SG24DTRATI0/TRATI0 

DCV=(0R1/C0SCDM)/F1'‘2“2*R(1)4DR1*F2/F4*1 •5*(1»1,5*F5/F4>)/Cy 
DRSQSDRJ*2*R<I) 

DDRSQ s-DRSQ* C DRHO+D VM 1 /VN< I ) > 


1930 
19406 
1950 
1963 
1970 
1980 
1990 
20006 
2010 
2020 
2030 
2040 
2030 10 

2060 
2070 
2075 
20 S0 
2090 
2100 


DDVMl=»F5*<DRH0vDRl/RCI )‘*DCV)-F6*DVHl/VM(n 
-LAM/Vtt(I >/R( 1 ) “2*< LAt4P2*nCU(I )*DRt /R (I ) ) 

VMlsVKCl l+CDVMH-DDVMI^DHFCI )/2>*DMF(l > 
Rl2(RCI)'2+(l)RSQ+DDRSQ*DMFa>/2)cDMFCI))'‘,5 
TRATIO=1-KV0*2>)‘RCUU4-1 I^ROTO-CRCUCl+l >/R 1 ) '^-VM 1 '‘2 >/C5 
RHOrRHO0*TRATIO'*G2*PLOS 

CVl = I /( 1 /CB+CRM-R I ) /CCS (8H) )+(R I *2-RM*2 )*SI MCBM) *2/ 
<Rr2+Rfr2*SI W{BH)''2) "I ,5 

DVM2= CV 1 /2/FI /RHO/R 1+LAM/VM !♦ (RCU < 1+ 1 ) /R I *2-C4 > 

DR2= .5/PI /BHO/VMl /R l*COS CBM ) 
VM(I+l)=VMa)4-(DVMI+DVM2)*DMFa)/2 
Ra+l)=R(I)+(0Rl+DR2)*0KFCI>/2 
CONTINUE 

* VMlTsVMCKS) J K=K+I 

V21=<CU0*2+VHn2)*2)'‘,9 
IF<X.GT.20) GO TO 32 
IF(ABSC1»V21/V2>.LT.. 00001)00 TO 18 
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APPENDIX III Ccont'd) 


Zlld 
2120 
2130 
2140 19 

2190 
2160 

2169 

2170 
2180 
2190 18 

2200 « 

2210 

2220 

2230 

2240 

2290 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340& 

2390 

2360 

2365 

2373 

2380 

2390& 

2400 

2410 

2420 

2430 

2440 

24906 

2460 

2470 

2480 

2490 

2900 20 

2910»t 
2920 
2930* 
2940 
2950 
2960 
2378 
2980 
2990 
2G00 
2610 22 
2620 
2630 


IFCK.NE.l) GU TO 19 
VM(KS)=,994ViKKS) 

60 TO 9 

V«<KS)=VWlT+<V«lT-V«2T)/(V21-V22>*(V2»V2i) 

DEC=VM<KS)/VH1T“1 

IF<ABS(DEC)oGT..2) VPKKS )sVNlT*( l+,2«:DEC/ABS (DEC ) > 
IF(KR,EQ,l)RM=RRt-<RK2-RMl>/(VW2T-VMm + (VMlT“VNU)) 

GO TO 9 

IF(KR.EQ.l) GO TO 21 

DO 20 J=l,6 
I=8-J 

TRATIO£l>KV0*2*RClKI)*ROTO- (RCUCI )/R(I ) > '“2-VMa ) '2)/C3 
RHO=RHO0»TRATI 0'‘G2«PL0S 
Fl = l/CM+<RM-Ra>>/COS{eR) 

F2SSIN(BW)'‘2 ; F5=Ra>^2-RK*2 
F4sR < I ) -‘2+ (R M4tSI N <BM ) ) *2 

CVsl/Fl+F2*F3/F4*1.9 t IF(ABS(CV) .LT.2.E- 18) CV=2.E-18 
F9=CV/2/PI/RH0/R(I) 

F6=LAM/VMa)*(RCUa)/R(I )"2*C4) 

DVWl s F5+F6 

0RU.9/PI /RK0/V«(I)/8a)*C0S(BM> 

DTRATI Os (- C4 #L AR4R CU a ) /R ( I > '^♦(LAM+RC U ( I) /R ( 1 )4DR 1 ) -VW< I ) VM I ) 
♦2/C3 

DRH0sQ2*DTRATl0/TRATI0 

DCVs (DRI /C03(BM) /FI '^2-2*R (I )*DR1»F2/F4'‘l .9*U“ 1 .9+F3/F4) )/CV 
DRSQsDRl*2>l‘Ra) 

DDftSQs»ORSQ*(DRHO+DVWl/V?Kl)) 
D0VMls-F9*(DRH0+DRl/Ra)-DCV)-F6*DVt1l/VH(I) 
-LAN/VM(I)/R(1)*2*(LAN+2*RCU(I)»DR1/Ra )) 

VIil = VMa)-(DVM 14 -DDVWl*DMF(l-l)/ 2 >*DI 1 F<I-l ) 
Rls(R(I)‘' 2 -*<DRSQ-DDRS 0 *DHF<I«l)/ 2 )»DHF(l-l )>*,9 
TRATIO 5 l+(V 0 ^ 2 *RCUCI+l)»ROTO-(RCUa+l)/Rl)^-VMl‘ 2 )/C 3 
R 1 {O=RHO 0 *TRATIO''G 2 *PLOS 

CW 1 s I /U /CM+<RM*R I ) /C OS (Bn ) >+ (R 1 *2-Rn *2 )*SI K (BPI ) *2/ 
(Rl'’2+Rn'*2*SlN(BM>'*2)*1.9 

OVM2sCVl/2/Pl/RHO/Rl+LAM/Vni*(RCU(I«l>/Rl'*2-C4) 

OR2s*9/PI/RHO/Vm/Rl*COS(Bf1> 

VMa-l)sVfl(l)<*(DVm+DVH2)*DMFa-l)/2 

Ra-l):R(l>-{DRl+DR2)4DMFa-l>/2 

CONTINUE 

IF(XR.EG*2) GO TO 21 

RIsR(12) J RKsRd) | PRlsV«(12) j PR2sVM(l> 

TRATIO=l+(V0“2*RCU(l)*ROTO-(RCU(l)/RU)) •‘2-VMU)'‘2)/C3 

PR5 sP 0*TRATIO*Q1«PLCS 

RK0lsRH0e«TRAn0'*'G2»PL0S 

PTINTS0 

DO 22 Is 1,11 

PnfJTsPTlBT+(PRATI 0 (l)+PRATI 0 a-H))/ 2 *(R(X+l)'“ 2 -R(I)* 2 ) 

CONTI WUE 

PTINTspnNT«PAf)*PLOS*Pl 
GO TO 59 
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APPENDIX III Ccont'd) 


2S40 

2650 21 DO 25 1=1,12 

2660 TRAT!O=l+<V0*2*RCUa)»ROTO-(RCUa)/R<I)>'‘2-VM<l)*2)/C3 

2670 P(I)=Pe#TRAT10‘*ai^PL0S 

26C0 I F(1 .EO. 1 >RHOl=RKO0*TRATIO'‘G2*PLOS 

2690 25 CQtlTU^UE 

2700 RHsRCl) I R«=RC7) | RTsR(l2) 

2710 DO 30 1=1,11 

2720 PRl=PRl+PI*(Ra+n-R(I))*(R(H-l )+R(I> )+(PAM-<P(n+P<I+| >)/2) 
,2730 PR3=PR5+Di1F<l)*(VMa)+VW<I+l))/2/U+R«*2*SIN(BM)''2/R<i>"‘2) 

27d0 30 CO»n»U£ 

2750 PR2=PI»!=RH'‘2+<PAf1-P<I>) 

2755 1F(ABS<KR-2KA8S<IHD-1)+ABS(R0T-1>)35.34,35 

2770 32 PRINT 33,KS,VM(KS),VW(12).R{KS),R<12) 

2775 QO TO 35 . 

2700 33 FORMATOH , R0TFL02 NOT CONVERGED K=20* ,11 ,2F7. 1 ,2F7.3) 

2785 34 PRINT 58 

27S0 DO 37 1=1,12 

2600 VI = VXU>«VIN 

2810 DTI = <1“EFF>*V0‘“2/C3*T0 

2815 TSI = T0-DII-’VI“2/C3>!<T0 

2820 PSlsPa*(TSX/<T0-DTl))'‘61 

2830 UB=Va*'2*R(I>/2/KCU0/EFF 

2840 TTR2=T0+(V0*2«RCUa)<i=ROTO/EFF-(RCUa>/R(I))*2+<UB-RCUa>/R(l>)*2 

284 1£ >/C34<T0 

2850 CPRs«Pa)/PSI>*(l/Gl>-l)/<TTR2/TSI-I) 

2852 SOL=SOLT+(SOLH-SOLT>*<R<12)/Ra)-1>/<R(12)/R<l)-1> 

2853 WI = <VI'^+UB'‘2)'‘.5 

2854 W0=(VB(i)'‘2+<HCU(I)/Ra))'‘2)‘*.5 

2855 DFACT=l-a0/WI+RCU(I)/2/S0L/R<I)/WI 

2860 PRINT 36,I,Ra),P<n-PAM,MI.W0,Vl.UB,S0L,CPR.DFACT 
2865 37 CONTINUE 

2870 36 FORMAKIH , I2,F8»3,F8*0,4F7.0,3F8.3) 

2875 38 FORMAT UH ) 

2880 35 RETURN {END 

2890* 

2900 SUBROUTINE INT2(X,Y,EI ,E2.«F,U 

2910 REALMF,t 

2920 DIMENSION X(3) ,Y(3>,E1 (3) ,E2(3) 

2930 DIMENSION C1(3),C2(3) 

2940 CALL COEF(X.Y,EI ,CI ) 

2950 CALL C0£RX.Y,E2,C2) 

2960 01=C1(2)*C2<3)-CH3)*C2(2) 

2970 D2=CI<5)*C2(I)-C1(I)*C2<3> 

2980 D3=C1(1}*C2<2)-C1(2)*C2(1) 

2990 DF1=1 t DF2=1 

3000 MFsX(n+DFl*(Dl/D3-X(l)) 

5010 L=Yn>+DF2*<D2/D3-Y<l>) 

3020 DEL1=1-MF/X<1) t DEL2=1-L/Y<1 ) 

3050 1F(ABS(DEH>.GT,.1> ttFs(| -.UABS (DELI) /DEL 1)*XU) 

3040 IF(ABS<0EL2).GT,.2) L=<l>.2*ABS<DEL2)/DEL2>*Y( 1> 

3050 RETURN f END 

3060* 


OF THE 
POOft 
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APPENDIX III (concluded) 


3070 SUBRODTIME C0EF(X,Y,E*C) 

32Se DlKEtiSIOH X<3) ,Y(3),S(3> 

30S0 DIMENSION C(3) 

3106 01=XU)+Y(2)+X<2)<=Y(3)+X(3)*YCI>»X<3)*Y<2>»X(2)*Y(I>-XU>*Y(3) 

31 10 D2=E<1)=^Y(2)+E(2>VY<5)+EC3)«YU )‘*E<3)*Y<2>»E(2>*YU )-E(l>*Y(3) 

3120 03:X<l)*E(2)'J-X(2>*EC3)+X(3)4<£(l)-X(3)*E(2)-X(2)*E<l)-Xn)*E(3) 

3150 D4=X(l)*Y€2)»E(3>+X(2)*Y(3)^Ea>+X<5)*Ya)*E(2) 

3I40& -XC3)*Y(2)«E<1>-X(2)*YU)*£C3>-X<1>*Y<3>*E(2> 

31&0 CC1)=D2/Dl 

3160 C(2)=D3/D1 

3170 CC3):D4/Di 

3180 HETUnti t END 

3190 FU«Cno» TAK(XXX>;7AU:Slti(XXX>/COS(XXX}$ RETURll } EKD 



TABLE I 



Outer Flowpath Inner Flowpath 

Flow* Exit Diameter Exit Diameter 


Model 

R/RT 

Angle 

In. 


In. 

iCml 

1 

0.4 

10 

3.116 

7.915 

1.246 

3.165 

2 

0.4 

0 

3.116 

7.915 

1.246 

3.165 

3 

0.4 

-10 

3.116 

7.915 

1,246 

3,165 

4 

0.4 

-20 

3.116 

7.915 

1.246 

3.165 

5 

0.55 

10 

3.420 

8,687 

1.879 

4,773 

6 

0.55 

0 

3.420 

8.687 

1.880 

4,773 

7 

0.55 

-10 

3.420 

8.687 

1.878 

4,773 

8 

0.55 

-20 

3.420 

8.687 

1.879 

4.773 

9 

0,70 

10 

4.000 

10.160 

2.795 

7.099 

10 

0.70 

0 

4.000 

10.160 

2.803 

7.099 

**E 

0.638 


3.710 

8.052 

2.366 

6.010 


* Radially out designated positive angle 
** Scaled from LF336/E lift fan (see Reference 2) 
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TABLE II 


Inner and Outer Exit Contour Geometries of the Base Pressure Models 


MODEL; 


2 


3 


4 


^^RT 









0 

0.400 

1.000 

0.400 

1.000 

0.400 

1,000 

0.400 

1.000 

- 0.128 

0.356 

0.983 

0.400 

1.000 

0.440 

1.016 

0.479 

1,033 

- 0.258 

0.306 

0.967 

0.400 

1.000 

0.478 

1.034 

0.549 

1.068 ’ 

- 0.385 

0.263 

0.951 

0.402 

1.003 

0.517 

1.053 

0.612 

1.105 

" 0*513 

0.237 

0.950 

0.411 

1.019 

0.563 

1,087 

0,674 

1.155 

- 0.642 

0.225 

0.962 

0.430 

1.048 

0.599 

1,137 

0.688 

1.206 

- 0.770 

0.229 

0.990 

0.459 

1.096 

0,625 

1.191 

0.703 

1.254 

- 0.899 

0.245 

1.033 

0.499 

1.155 

0.639 

1.239 

0,703 

1.284 

- 1.027 

0.278 

1.091 

0.548 

1.205 

0.642 

1.284 

0.689 

1.284 

- 1.155 

0.327 

1.151 

0.589 

1.253 

0.642 

1,284 

0.664 

1.284 

- 1.284 

0.387 

1.202 

0.618 

1.284 

0.642 

1,284 

0,648 

1.284 

- 1.412 

0.447 

1.249 

0.635 

1.284 

0.642 

1.284 

0.642 

1.284 

- 1.540 

0.507 

1.284 

0.642 

1.284 

0.642 

1.284 

0,642 

1.284 

- 1.669 

0.558 

1.284 

0.642 

1.284 

0.642 

1.284 

0.642 

1.284 

- 1.797 

0.596 

1.284 

0.642 

1.284 

0.642 

1,284 

0.642 

1.284 

- 1.926 

0.623 

1.284 

0.642 

1.284 

0.642 

1,284 

0.642 

1.284 

- 2 . 054 

0.637 

1.284 

0.642 

1.284 

0,642 

1,284 

0,642 

1.284 

- 2.182 

0.642 

1.284 

0.642 

1.284 

0,642 

1.284 

0.642 

1,284 



TABLE II (continued) 


UODEL: 5 g 


Z/jjT 


(R/rt>t 

(R/rt>h 


<«/rt>h 

(R/rt>t 


RT T 

0 

0.550 

1.000 

0.550 

1.000 

0.550 

1.000 

0.550 

1,000 

-0.X28 

0.520 

0.984 

0.550 

1.000 

0.580 

1.016 

0,608 

1,033 

-0.258 

0.488 

0.968 

0.550 

1.000 

0.608 

1.034 

0,662 

1,065 

-0.385 

0.457 

0.953 

0.551 

1.001 

0.635 

1.053 

0,712 

1.101 

-0.513 

0.435 

0.950 

0.560 

1.015 

0.649 

1.085 

0.748 

1.145 

-0.642 

0.422 

0.961 

0.575 

1.042 

0.651 

1.127 

0.769 

1.170 

-0.770 

0.418 

0.984 

0.584 

1.083 

0.639 

1.170 

0.777 

1,170 

-0,899 

0.422 

1.023 

0.585 

1.126 

0.616 

1.170 

0.771 

1.170 

-1 . 027 

0.436 

1.073 

0.585 

1.169 

0.596 

1.170 

0,753 

1,170 

-1.155 

0.458 

1.117 

0.585 

1.170 

0.585 

1.170 

0,720 

1.170 

-1,284 

0.490 

1.160 

0.585 

1.170 

0,585 

1.170 

0,674 

1.170 

-1,412 

0.528 

1.170 

0,585 

1.170 

0.585 

1.170 

0.636 

1.170 

-1,540 

0.557 

1.170 

0.585 

1.170 

0.585 

1.170 

0,609 

1,170 

-1,669 

0.575 

1.170 

0,585 

1.170 

0.585 

1.170 

0.592 

1.170 

-1,797 

0,584 

1.170 

0.585 

1.170 

0.585 

1.170 

0.585 

1.170 

-1.926 

0,585 

1.170 

0.585 

1.170 

0.585 

1.170 

0.585 

1.170 

-2 . 054 

0.585 

1.170 

0.585 

1.170 

0.585 

1.170 

0,585 

1,170 

-2.182 

0.585 

1.170 

0.585 

1,170 

0,585 

1.170 

0.585 

1.170 



MODEL : 


9 




0 

0.700 

1.000 

- 0.128 

0.679 

0.985 

- 0.258 

0.657 

0.970 

- 0,385 

0.633 

0.957 

- 0.513 

0.598 

0.953 

- 0.642 

0.560 

0.961 

- 0.770 

0.532 

0.981 

- 0.899 

0.513 

1.000 

- 1.027 

0.502 

1.000 

- 1.155 

0.500 

1.000 

- 1.284 

0.500 

1.000 

- 1.412 

0.500 

1.000 

- 1.540 

0.500 

1.000 

- 1.669 

0.500 

1.000 

- 1.797 

0.500 

1.000 

- 1.926 

0.500 

1.000 

- 2.054 

0.500 

1.000 

- 2.182 

0.500 

1.000 


ILE II ( continued ) 


10 E 






0.700 

1.000 

0.638 

1.000 

0.700 

1,000 

0.595 

1.029 

0.700 

1.000 

0.564 

1.058 

0.699 

1.000 

0.545 

1.078 

0,688 

1.000 

0.539 

1.078 

0.665 

1.000 

0.539 

1.078 

0.629 

1.000 

0.539 

1.078 

0.585 

1.000 

0.539 

1,078 

0.550 

1.000 

0.539 

1.078 

0.525 

1.000 

0.539 

1.078 

0.509 

1.000 

0.539 

1.078 

0,501 

1.000 

0.539 

1.078 

0.500 

1.000 

0.539 

1.078 

0.500 

1.000 

0.539 

1.078 

0.500 

1.000 

0.539 

1 .078 

0.500 

1,000 

0.539 

1.078 

0.500 

1.000 

0.539 

1.078 

0.500 

1.000 

0.539 

1.078 



U1 

■j4 


PRESSURE RATIO 1.3 

TRAVERSE AXIAL (IN.) 0.1 

DISTANCE (CM.) 0.25A 

BLOCKAGE 

MODEL: A1 X 

A2 X 

A3 X 

A4 X 

A5 X 

A6 X 

A7 X 

A8 X 

A9 X 

AlO X 

AE X 

51 X 

52 X 

53 X 

54 X 

55 X 

56 X 

57 X 

58 X 

59 X 

SIO X 

SE X 

SA2 
SA6 

SAIO X 


X X 


TABLE III 




TABLE IV 


Performance Comparisons for Axial ^ Swirling 
and Skev?ed Flow Distributions 


Average 

Base 

Pressure 

Coefficient 

Model c 

P2 


A2 

-.166 

S2 

-.840 

SA2 

-.484 

A6 

-.195 

S6 

-.740 

SA6 

-.341 

AIO 

-.211 

SIO 

-.500 

SAIO 

-.240 


Thrust 

Flow 

Coefficient 

Coefficient 

Crr 

c 

T 

F 

.979 

.972 

.695 

.827 

.903 

.954 

.964 

,991 

.684 

.888 

.981 

.993 

.888 

.980 

.660 

.942 

.861 

1.009 
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TABLE V 


Comparison of Test and Theory for Configurations 



SA2, SA6, SAIO. 

SE 




Average 

Base 

Pressure 

Thrust 

Flow 



Coefficient 

Coefficient 

Coefficient 

Model 





SA2 

-.484 

.903 

.954 

Test 


-.463 

.910 

1.006 

Theory 

SA6 

-.341 

.981 

.993 

Test 


-.400 

.871 

1.008 

Theory 

SAIO 

-.240 

.861 

1.009 

Test 


-.272 

.858 

1.019 

Theory 

SE 

-.535 

.688 

.868 

Test 


-.553 

.628 

.941 

Theory 
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Figure 1. Analytical Model 


Case 1 , 


Parallel Streamlines and rm^O 



" (R - Rm)/cos 


Case 2 . Constant Area Between Streamlines and rjjj = 0 



Figure 2. Streamline Curvature Approximations 
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Figure 5, Photograph of the Base Pressure Model Assembly 


1 OF THE 
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Figure 6, Photograph of the Axial and Swirl Vane Assemblies 
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Inner Fl owpath 
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Figure 7, 


Hub Base Pressure Test Model Schematic 



Figure 8. 


Photograph of 


the Base Pressure Model Installed in 
the Test Stand 
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45° 


0.305) 


Figure 9. Schematic of Exhaust Traverse Probe 







Pressure Coefficient (CPT^ CPS) 




Figure H. Total and Static Pressure Coefficient 

Distribution for Models 1 and 2, Axial Flow 
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Pressure Coefficient (CPT, CPS) 




Figure 12. Total and Static Pressure Coefficient 

Distribution for Models 3 and 4, Axial Flow 
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Pressure Coefficient (CPT, CPS) 




Figure 14. Total and Static Pressure Coefficient 

Distribution for Models 7 and 8, Axial Flow 
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Pressure Coefficient (CPT, CPS 



Radius (R/R,p) 

Figure 15. Total and Static Pressure Coel'iicieiit 

Distribution for Models 9 and 10, Axial Flow 
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Pressure Coefficient (CPT, CPS) 



Figure 16. Total and Static Pressure Coefficient 
Distribution for Model E, Axial Flow 
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Radius (R/R.J,) 

Figure 17, Total and Static Pressure Coefficient 

Distribution for Models 1 and 2, Axial Flow 
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Radius (R/R^) 

Figure 18. Total and Static Pressure Coefficient 

Distribution for Models 3 and 4, Axial Flow 





Pressure Coefficient (CPT, CPS) 




Figure 19. Total and Static Pressure Coefficient 

Distribution for Models 5 and 6, Axial Flow 
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Radius (R/El^) 

Figure 20. Total and Static Pressure Coefficient 

Distribution for Models 7 and 8, Axial Flow 
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Pressure Coefficient (CPT, CPS) 




Figure 21. 


Total and Static Pressure Coefficient 
Distribution for Models 9 and 10, Axial Flow 







Pressure Coefficient (CPT, CPS) 



Figure 22. Total and Static Pressure Coefficient 
Distribution for Model E, Axial Flow 
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Axial Distance (Z/R^) 


Figure 29 . Flow Entrainment Chara 


Flow Ratio 



for Models 1 through 4, Axial Flow 





Flow Ratio 



Figure 30. Flow Entrainment Characteristics for Models 5 through 8, Axial Flow 
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Figure 32 , Wake Mixing Loss Characterist 


Wake Mixing Loss Coefficient (o)) 






Wake Mixing Loss Coefficient (cd) 
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Average Hub Base Pressure Coefficient (CPB) 


+ 0,2 



Flow Angle ( 0 ), Degrees 

Figure 35, Effects of Geometry on Average Hub Base 
Pressure^ Axial Flow 


-9 










Figure 38, Effects of Geometry on Flow Coefficient, 
Axial Flow 
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Pressure Coefficient (CPT, CPS) 



Radius (R/R^) 

Figure 44 Total and Static Pressure Coefficient 

and Swirl Angle Distribution for Models 
E, Swirling Flow, Close-Up Traverse 
Position 
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Figure 45. Total and Static Pressure Coefficient and 

Swirl Angle Distribution for Models 1 and 2, 
Swirling Flow, Downstream Traverse Location 
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Figure 46. Total and Static 
Swirl Angle Dist 


Swirl Angle (O'), Degrees 
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Figure 47. Total and Static Pressure Coefficient and 

Swirl Angle Distribution for Models 5 and 6, 
Swirling Flow, Downstream Traverse Location 
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Swirl Angle (o'). Degrees 




Pressure Coefficient (CPT, CPS) 



Figure 48. Total and Static Pressure Coefficient and 

Swirl Angle Distribution for Models 7 and 8, 
Swirling Flow^ Downstream Traverse Location 
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Swirling Angle (a). Degrees 





Pressure Coefficient (CPT, CPS) 



Figure 49. Total and Static Pressure Coefficient and 

Swirl Angle Distribution for Models 9 and 10, 
Swirling Flow, Downstream Traverse Location 
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Swirling Angle (O'), Degrees 




Pressure Coefficient (CPT, CPS) 


240 



Figure 50. Total and Static Pressure Coefficient and 
Swirl Angle Distribution for Model E, 
Swirling Flow, Downstream Traverse Location 
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Figure 55. Graphical Representation of the Flow Field for 
Models 9 and 10, Swirling Flow 
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Figure 58. Flow Entrainment Characteristics for Models 1 through 4, Swirling Flow 
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Average Hub Base Pressure Coefficient (C 
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Figure 63 . Effects of Geometry on Average Hub Base 
Pressure, Swirling Flow 
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Hub Base Pressure Coefficient (CP 
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Figure 64. Effects of Geometry on the Hub Outer Radius Base 
Pressure Coefficient, = 10 & 0, Swirling Flow 
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Hub Base Pressure Coefficient (CPj^ 



Axial Distance (Z/R^) 


Figure 65. Effects of Geometry on the Hub Outer Radius Base 

Pressure Coefficient, ^ = -10 & -20, Swirling Flow 
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Main Flow Thrust Coefficient (c 



Figure 66. Effects of Geometry on Main Flow Thrust, 
Swirling Flow 
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Flow Angle (B), Degrees 

Figure 68. Effects of Geometry on Flow Coefficient 
Swirling Flow 
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Figure 70* Effects of Pressure Ratio on Total and Static 
Pressure Coefficients and Swirl Distribution 
for Model 6^ Swirling Flow 







Figure 71, Effects of Pressure Ratio on Total and Static 
Pressure Coefficients and Swirl Distribution 
for Model 10, Swirling Flow 
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Figrure 72 , Effects of Pressure Ratio on the Average Hub 
Base Pressure Coefficient, Swirling Flow 


- 131 - 



Hub Base Pressure Coefficient (C 




1.0 



0.9 




0.7 

1.2 

Pressure 

Figure 74. Effects of Pres 
Thrust Coefficf 
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Figure 77 • Photograph of the Blockage System Installed in the 

Test Stand 
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Figure 78. Effects of Blockage on Total and Static Pres- 
sure Coefficients for Model 2^ Axial Flow 
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Figure 79 ♦ Effects of Blockage on Total and Static Pressure 
Coefficients for Model 6^ Axial Flow 
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Figure 80. Effects of Blockage on Total and Static Pres- 
sure Coefficients for Model 10, Axial Flow 
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Figure 82. Effects of Blockage on Total and Static Pres- 
sure Coefficients and Swirl Distribution for 
Model 2, Swirling Flow 
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Figure 86, Effects of Blockage on The Flow Coefficient, 
Swirling Flow 
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Figure 88, Total and Static Pressure Coefficient 

and Swirl Angle Distribution for Model 6, 
Skewed Flow Distribution 
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Figure 90. Angular Momentum Distribution From Test 
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Figure 91, W^ake Area Coefficient 
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Flow Angle (P)^ Degrees 

Figure 94. Comparison of Theoretical Predictions to Test Results 

for the Base Pressure^ Flow and Thrust Coefficients^ 

R /R =0.55^ Swirling Flow 
11 1 
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Figure 95, 


Comparison of Theoretical Predictions to Test Results 
for the Base Pressure, Flow and Thrust Coefficients, 

= 0.70^ Swirling Flow 
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Figure 96 , Comparison of Theoretical Predictions to Test Results 
for the Base Pressure^ Flow and Thrust Coefficients, 
R^/R^ = 0.40, Axial Flow 
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Figure 97. Comparison of Theoretical Predictions to Test Results 
for Base Pressure^ Flow and Thrust Coefficients, 

Rjj/R^ = 0.55, Axial Flow 
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Figure 9a Comparison of Theoretical Predictions to Test Results 
for Base Pressure, Plow and Thrust Coefficients, 

Rjj/R^ = 0.70, Axial Flow 
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